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ACTIVE  AND  PASSIVE  REMOTE  SENSING  OF  ICE 


ANNUAL  REPORT 


This  annuikl  report  covers  research  under  the  sponsorship  of  the  ONR  contract 
N00014-89-J-1107  from  October  1,  1990  to  September  30,  1991.  We  have  published  15 
journal  and  conference  papers  and  3  student  theses  during  this  period. 

1.  Multi-layer  Random  Medium  Model  for  Sea  Ice 

✓ 

N  ■ 

Fully  polarimetric  scattering  of  electromagnetic  waves  from  snow  and  ice  is  stud¬ 
ied  with  a  multi-layered  random  medium  model  and  applied  to  interpret  experimental 
data  obtjiined  under  laboratory  controlled  conditions  such  as  CRRELEX.  The  snow  layer 
is  modeled  as  an  isotropic  random  medium.  The  sea  ice  is  described  as  an  anisotropic 
random  medium  due  to  the  nonspherical  shape  of  brine  inclusions.  The  underlying  sea 
water  is  considered  as  a  homogeneous  half-space.  The  random  media  in  both  layers  are 
characterized  by  three-dimensional  correlation  functions  with  varieinces  and  correlation 
lengths  corresponding  to  the  fluctuation  strengths  and  the  physical  geometries  of  the  inho¬ 
mogeneities,  respectively.  The  strong  fluctuation  theory  is  used  to  calculate  the  effective 
permittivities  of  the  random  media.  The  distorted  Born  approximation  is  then  employed 
to  obtain  the  covariance  matrix  which  represents  the  fully  polarimetric  scattering  prop¬ 
erties  of  the  snow-ice  media.  It  has  been  shown  that  the  polarimetric  covariance  matrix 
contains  more  information  than  the  conventional  scattering  coefficients  on  the  remotely 
sensed  media. 

In  saline  ice  sheets  under  quiescent  condition,  the  background  ice  grows  in  colum¬ 
nar  form  and  saline  water  is  trapped  between  ice  platelets  in  the  form  of  brine  inclusions 
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which  are  usually  ellipsoidal.  The  ice  tends  to  grow  vertically  downward  rendering  the 
ellipsoidal  inclusions  aligned  preferably  in  the  vertical  direction  and  the  crystallographic 
C  axes  parallel  to  the  horizontal  plane.  In  this  case,  the  C  axes  are,  however,  random  in 
azimuthal  direction.  The  strong  fluctuation  theory  is  extended  to  account  for  vertically 
aligned  ellipsoideil  brine  inclusions  with  C  axes  randomly  oriented  in  the  horizont^ll  direc¬ 
tion.  The  brine  inclusions  are  described  by  three-dimensional  local  correlation  functions. 
The  configuration  average  over  the  azimuthal  orientation  angles  is  carried  out  in  the  pro¬ 
cess  of  deriving  the  global  correlation  tensor.  The  distorted  Born  approximation  is  applied 
to  obtain  the  covariance  matrix  for  the  multi-layered  snow-ice  configuration.  The  theo¬ 
retical  results  show  non-zero  cross-polarized  returns  under  the  first-order  distorted  Born 
approximation.  We  have  also  compared  the  results  with  experimental  data  obtained  by 
the  US  Army  Cold  Regions  Research  and  Engineering  Laboratory  (CRREL). 

In  sea  ice,  the  scatterers  can  have  various  shapes,  sizes,  and  permittivities.  We 
have  also  investigated  the  modeling  of  radar  backscatter  from  random  media  with  multiple 
scatterer  species.  We  consider  each  type  of  scatterers  as  a  species  which  can  take  on  a  shape, 
size,  and  complex  permittivity  different  from  other  species.  The  multiple  species  in  the 
random  medium  are  considered  as  randomly  oriented  ellipsoids  and  described  by  multiple 
three-dimensional  ellipsoidal  local  correlation  functions.  The  variances  and  correlation 
lengths  of  the  correlation  functions  characterize  the  fluctuation  strengths  and  the  physical 
geometries  of  each  species  of  scatterers.  The  effective  permittivity  of  the  random  medium 
is  derived  under  the  strong  fluctuation  theory  and  the  polarimetric  scattering  coefficients 
are  calculated  for  the  layer  configuration  with  the  distorted  Born  approximation.  Due 
to  the  non-spherical  shape  and  the  random  orientation  of  the  scatterers,  the  correlation 
coefficient  between  the  HH  and  W  returns  has  a  magnitude  different  from  unity  and  a 
small  phase  angle.  The  scattering  coefficients  are  also  used  to  calculate  the  Mueller  matrix 
for  synthesis  of  polarization  signatures.  The  co-polarized  signature  of  the  random  medium 
has  a  rather  straight  distortion  track  and  a  recognizable  pedestal. 
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2.  Polarimetric  Radar  Calibration 


Accurate  calibration  of  polarimetric  radar  systems  is  essential  for  the  polarimetric 
remote  sensing  of  earth  terrain.  Polarimetric  calibration  algorithms  using  in-scene  re¬ 
flectors  and  distributed  t2U'gets  are  studied.  The  transmitting  and  receiving  ports  of  the 
polarimetric  radar  are  modeled  by  two  unknown  polarization  transfer  matrices.  These 
unknown  matrices  are  determined  using  the  the  measured  scattering  matrices  from  the 
calibration  targets.  For  the  case  of  polarimetric  calibration  using  three  in-scene  reflec¬ 
tors,  a  Polarization-basis  Transformation  (PT)  technique  is  introduced  to  find  out  a  new 
transmitting  and  receiving  polarization  basis  under  which  the  scattering  matrices  of  the 
calibration  targets  will  fall  into  one  of  six  simpler  sets.  The  calibration  solution  can  then 
be  solved  easily  in  the  new  polarization  basis  and  converted  to  obtain  the  solution  in  the 
original  polarization  basis.  The  uniqueness  of  polarimetric  calibration  using  three  targets 
is  addressed  for  aU  possible  target  combinations.  The  PT  technique  can  also  be  appUed 
to  the  polarimetric  calibration  using  a  combination  of  arbitreiry  in-scene  reflectors  amd 
distributed  targets.  The  effect  of  misalignment  of  calibration  targets  and  the  sensitivity 
of  polarimetric  calibration  algorithms  to  the  noise  are  illustrated  by  investigating  several 
sets  of  Ccilibration  targets. 

3.  Classification  Studies  of  Polarimetric  SAR  Images 

We  have  studied  the  SAR  image  classification  by  using  the  neural  network  meth¬ 
ods.  Supervised  methods,  including  both  conventional  Maximum  Likelihood  (ML)  2ind 
more  recent  multi-layer  perceptron  neural  network  classifiers  have  yielded  higher  accuracy 
than  unsupervised  techniques,  but  suffer  from  the  need  for  human  interaction  to  prede¬ 
termine  classes  and  training  regions.  In  contrast,  unsupervised  methods  determine  classes 
automatically,  but  generally  show  limited  ability  to  accurately  divide  terrain  into  natural 
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cleisses.  We  introduced  a  new  terrain  classification  technique  to  discriminate  sea  ice  sig¬ 
natures  in  polarimetric  SAR  images  by  utilizing  unsupervised  neural  networks  to  provide 
automatic  classification,  but  employing  an  iterative  algorithm  which  overcomes  the  poor 
accur2u:y  of  other  unsupervised  techniques. 

Several  types  of  unsupervised  neural  networks  are  first  applied  to  the  cl2issification 
of  SAR  images,  and  the  resiilts  are  compared  with  those  of  more  conventional  unsupervised 
methods.  Neural  network  approaches  include  Adaptive  Resonance  theory  (ART),  Learning 
Vector  Quamtization  (LVQ),  and  Kohonen’s  self-organizing  feature  map.  Conventional 
classifiers  utilized  aire  the  migrating  means  clustering  algorithm  aind  the  K-means  clustering 
method.  Preprocessing  is  performed  with  the  SAR  images  to  reduce  speckle  noise  and 
stabilize  the  training  process  for  the  neural  networks.  Results  after  preprocessing  show  that 
LVQ  and  Kohonen’s  self-organizing  feature  map  outperform  the  conventional  unsupervised 
classifiers,  but  are  still  inferior  to  supervised  methods. 

To  overcome  this  poor  accuracy,  an  iterative  algorithm  is  constructed  where  the 
SAR  image  is  reclassified  using  a  Maximum  Likelihood  (ML)  classifier.  Training  of  the  ML 
classifier  is  performed  using  a  training  data  set  first  classified  by  the  above  unsupervised 
method,  thus,  requiring  no  humsui  intervention,  and  preserving  the  unsupervised  nature 
of  the  overciU  classification  scheme.  The  process  is  then  repeated  iteratively,  training  a 
second  ML  classifier  using  data  classified  by  the  first.  It  is  shown  that  this  algorithm 
converges  rapidly,  and  significantly  improves  classification  accuracy.  Performance  after 
convergence  is  seen  to  be  comparable  to  that  obtained  with  a  supervised  ML  classifier, 
while  maintaining  the  advantages  of  an  unsupervised  technique. 

The  unsupervised  and  iterative  techniques  developed  have  been  applied  to  the 
polarimetric  SAR  images  of  Beaufort  sea  ice  acquired  by  the  C-,  L-,  and  P-band  SAR 
instruments  of  Jet  Propulsion  Laboratory.  The  results  obtained  with  the  new  algorithms 
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Me  compMed  with  the  results  obtained  with  other  techniques  by  classifying  terrain  features 


in  polarimetric  SAR  images. 


4.  K-Dbtribution  and  Polarimetric  Terrain  Radar  Clutter 

A  multivariate  K- distribution  has  been  developed  to  model  the  statistics  of  fully 
polarimetric  radar  data  from  earth  terrain  with  polarizations  HH,  HV,  VH,  and  W.  In 
this  approach,  correlated  polarizations  of  radar  signals,  as  characterized  by  a  covariance 
matrix,  are  treated  as  the  sum  of  N  n-dimensional  random  vectors;  N  obeys  the  negative 
binomial  distribution  with  a  parameter  a  and  mean  N.  Subsequently,  an  n-dimensional 
K-distribution,  with  either  zero  or  nonzero  mean,  is  developed  in  the  limit  of  infinite  N  or 
illuminated  area.  The  probability  density  function  (PDF)  of  the  K-distributed  vector  nor¬ 
malized  by  its  Euclidean  norm  is  independent  of  the  parameter  a  and  is  the  same  as  that 
derived  from  a  zero-mean  Gaussian-distributed  random  vector.  The  above  model  is  well 
supported  by  experimental  data  provided  by  MIT  Lincoln  Laboratory  and  the  Jet  Propul¬ 
sion  Laboratory  in  the  form  of  polarimetric  measurements.  The  results  are  illustrated  by 
comparing  the  higher-order  normalized  intensity  moments  and  cumulative  density  func¬ 
tions  (CDF)  of  the  experimental  data  with  theoretical  residts  of  the  K-distribution. 

5.  Composite  Volume  and  Surface  Scattering  Model 

Among  the  various  theoretical  models  applied  to  study  the  electromagnetic  wave 
scatterings  from  geophysical  terrain,  such  as  snow  and  ice,  the  radiative  transfer  theory 
has  drawn  intensive  attention  in  the  microwave  remote  sensing  society  during  the  past 
years.  In  most  of  the  scattering  models,  the  volume  scattering  and  the  surface  scattering 
effects  have  been  investigated  separately.  Recently,  there  has  been  a  growing  interest  in 
the  construction  of  composite  models  which  can  take  into  account  both  types  of  scattering. 


We  derived  the  first  order  iterative  solution  to  the  vector  radiative  transfer  equations  for  a 
two-layer  medium  with  a  diffuse  top  boundary  and  an  irregular  bottom  boundary  of  Gaus- 
siein  roughness.  The  Kirchhoff  approximation  and  the  geometrical  optics  approach  with 
shadowing  correction  are  used  in  formulating  the  boundary  conditions.  To  demonstrate 
the  utilities  of  the  theory,  randomly  oriented  spheroidal  discrete  scatterer  model  is  used  to 
calculate  the  backscattering  coefficients  from  soybean  field  in  different  growing  stages  and 
compared  to  the  experimental  measurements.  Good  agreement  has  been  achieved  for  both 
the  co-polarized  and  the  cross-polarized  data.  It  is  observed  that  the  presence  of  the  rough 
surface  can  significantly  enhance  the  backscattering  at  small  incident  angles  and  increase 
the  cross-polarized  returns.  The  polarization  signatures  calculated  based  on  the  Mueller 
matrix  show  a  straight  distortion  track  and  an  observable  pedestal.  Numerical  comparison 
to  the  backscattering  coefficients  calculated  by  using  planar  bottom  boundary  conditions 
with  or  without  the  incoherent  addition  of  the  rough  surface  effects  are  iJso  made. 

6.  Polarimetric  Passive  Remote  Sensing 

The  concept  of  polarimetry  in  active  remote  sensing  is  extended  to  passive  remote 
sensing.  The  potentizd  use  of  the  third  and  fourth  Stokes  parameters  U  and  V ,  which 
play  an  important  role  in  polarimetric  active  remote  sensing,  is  demonstrated  for  passive 
remote  sensing.  It  is  shown  that,  by  the  use  of  the  reciprocity  principle,  the  polarimetric 
parameters  of  passive  remote  sensing  can  be  obtained  through  the  solution  of  the  associated 
direct  scattering  problem.  In  particular,  the  full  polarimetric  information,  including  the 
corresponding  brightness  temperatures  of  U  and  V,  can  be  obtained  from  the  solution  of 
the  direct  scattering  problem  for  four  different  polarizations  of  the  incident  wave.  These 
ideas  are  applied  to  study  polarimetric  passive  remote  sensing  of  periodic  surfaces.  The 
solution  of  the  direct  scattering  problem  is  obtained  by  an  integrid  equation  formulation. 
Incidence  on  a  penetrable,  lossy,  medium  is  considered.  Since  the  kernels  of  the  integral 
equations  are  the  periodic  Green’s  functions  and  their  normal  derivatives  on  the  surface. 
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rapid  evaluation  of  the  slowly  convergent  series  associated  with  these  functions  is  observed 
to  be  critical  for  the  feasibility  of  the  method.  The  study  has  shown  that  the  brightness 
temperature  of  the  Stokes  parameter  U  can  be  significant  in  passive  remote  sensing.  Values 
as  high  as  50  K  are  observed  for  certain  configurations. 

To  demonstrate  the  use  of  polarimetry  in  passive  remote  sensing  of  azimuthally 
asymmetric  features  on  a  terrain  surface,  an  experiment  was  designed  and  implemented.  A 
triangulax  corrugation  pattern  was  made  on  the  sandy  soil  surface.  Polarimetric  brightness 
temperatures  are  measured  with  horizontal,  vertic<il,  and  45°  polarization  orientations  for 
vzirious  observation  angles.  From  the  measured  temperatures,  absolute  values  as  high  as 
30— 40iir  of  the  third  Stokes  brightness  temperatures  are  observed.  A  theoretical  analysis 
of  the  data  indicates  that  the  high  values  of  XJ  are  caused  by  the  azimuthal  asymmetry 
on  the  remotely  sensed  soil  surface.  It  is  also  observed  from  the  experiment  that  the 
brightness  temperatures  for  all  three  Stokes  parameters  vary  as  the  observation  direction 
varies  from  being  parallel  to  the  surface  row  structure  to  being  perpendiculzir  to  the  row 
structure.  The  significant  implication  of  this  experiment  is  that  the  surface  asymmetry 
czm  be  detected  with  a  measurement  oiU  at  a  single  azimuthal  angle. 
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Polarimetric  Passive  Remote  Sensing  of  Periodic  Surfaces 
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Cambridge,  MA  02139,  USA 

Abstract-  The  concept  of  polarimetry  in  active  remote  sensing  is  extended  to  passive 
remote  sensing.  The  potential  use  of  the  third  and  fourth  Stokes  parameters  U  and 
V,  which  play  an  important  role  in  polarimetric  active  remote  sensing,  is  demonstrated 
for  passive  remote  sensing.  It  is  shown  that,  by  the  use  of  the  reciprocity  principle,  the 
polarimetric  parameters  of  passive  remote  sensing  can  be  obtained  through  the  solution 
of  the  associated  direct  scattering  problem.  These  ideas  are  applied  to  study  polarimetric 
passive  remote  sensing  of  periodic  surfaces.  The  solution  of  the  direct  scattering  problem 
is  obtained  by  an  integral  equation  formulation  which  involves  evaluation  of  periodic 
Green’s  functions  and  normal  derivative  of  those  on  the  surface.  Rapid  evaluation  of  the 
slowly  convergent  series  associated  with  these  functions  is  observed  to  be  critical  for  the 
feasibility  of  the  method.  New  formulas,  which  are  rapidly  convergent,  are  derived  for  the 
calculation  of  these  series.  The  study  has  shown  that  the  brightness  temperature  of  the 
Stokes  parameter  U  can  be  significant  in  passive  remote  sensing.  Values  as  high  as  60  K 
are  observed  for  certain  configurations. 

L  INTRODUCTION 

Although  the  tue  of  polarimetry  is  well  understood  in  active  remote  sensing  [1-4], 
little  attention  has  been  given  to  it  in  passive  remote  sensing.  In  active  remote 
sensing,  additional  polarimetric  information  is  obtained  from  the  third  and  fourth 
Stokes  parameters  U  and  V.  They  are  defined  as  2Re{EvE^)  and  2lm{EvE^) 
respectively,  where  the  brackets  denote  ensemble  averaging.  Therefore,  U  and 
V  are  the  measure  of  correlation  between  the  horizontal  and  vertical  components 
{Ej^  and  Ev)  of  the  electric  field  and  thereby  carry  further  information.  Tsang  [5] 
extended  the  concept  of  polarimetry  to  passive  remote  sensing  and  considered 
scattering  by  discrete  scatters  and  rough  surfaces  when  the  thermal  emission  of 
third  and  fourth  Stokes  parameters  are  nonzero.  Passive  polarimetric  measure¬ 
ments  have  been  extensively  used  in  radio  astronomy  to  study  planets  and  stars 
|6].  In  this  paper,  we  discuss  polarimetric  passive  remote  sensing  of  periodic  rough 
surfaces.  In  particular,  we  calculate  the  brightness  temperatures  of  U  and  V  as 
can  be  obtained  by  passive  measurements. 

The  reciprocity  principle  is  used  to  relate  active  remote  sensing  parameters  to 
polarimetric  passive  remote  sensing.  The  emissivity  of  the  surface  for  a  particular 
polarization  can  be  obtained  as  one  minus  the  reflectivity  [7,8].  In  order  to  be  able 
to  get  the  full  polarimetric  information,  including  the  corresponding  brightness 
temperatures  of  U  and  V,  we  need  the  emissivities  of  four  different  polarizations: 
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To  calculate  the  scattered  fields,  the  y-component  of  the  surface  electric  and  P 
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all  the  remaining  components  of  the  field  can  be  determined  through  formulas  like  *  '  of  the  problem. 
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Classification  of  terrain  cover  using  polarimetric  radar  is  an  area  of  considerable  cur¬ 
rent  interest  and  research.  A  number  of  methods  have  been  developed  to  classify  ground 
terrain  types  from  fully  polarimetric  synthetic  aperture  radar  (SAR)  images,  and  these 
techniques  are  often  grouped  into  supervised  and  unsupervised  approaches.  Supervised 
methods,  including  both  conventional  Maximum  Likelihood  (ML)  and  more  recent  Multi¬ 
layer  Perceptron  classifiers,  have  yielded  higher  accuracy  than  unsupervised  techniques, 
but  suffer  from  the  need  for  human  interaction  to  predetermine  classes  and  training  re¬ 
gions.  In  contrast,  unsupervised  methods  determine  classes  automatically,  but  generally 
show  limited  ability  to  accurately  divide  terrain  into  naturcd  classes.  In  this  paper,  a 
new  terrain  classification  technique  is  introduced,  utilizing  unsupervised  neural  networks 
to  provide  automatic  classification,  but  employing  an  iterative  algorithm  which  overcomes 
the  poor  accuracy  of  other  unsupervised  techniques. 

Several  types  of  unsupervised  neural  networks  are  first  applied  to  the  classification  of 
SAR  images,  and  the  results  are  compared  with  those  of  more  conventional  unsupervised 
methods.  Neural  Network  approaches  include  Adaptive  Resonance  theory  (ART),  Learning 
Vector  Quantization  (LVQ),  and  Kohonen’s  self-organizing  feature  map.  Conventional 
classifiers  utilized  are  the  migrating  means  clustering  algorithm  and  the  K-means  clustering 
method.  With  both  neural  network  and  convention^  classifiers,  preprocessing  is  performed 
to  reduce  speckle  noise  and  to  stabilize  the  training  process.  Results  show  that  LVQ  is 
the  best  of  the  neural  network  techniques,  and  that  this  method  outperforms  all  of  the 
conventional  unsupervised  classifiers.  The  accuracy  of  even  the  LVQ  technique,  however, 
is  seen  to  remain  below  that  of  supervised  methods. 

To  overcome  this  poor  accuracy,  an  iterative  algorithm  is  proposed  in  which  the  SAR 
image  is  reclassified  using  a  Maximum  Likelihood  (ML)  classifier.  Training  of  the  ML 
classifier  is  performed  using  a  training  data  set  first  classified  by  the  above  unsupervised 
method,  thus,  requiring  no  human  intervention,  and  preserving  the  unsupervised  nature  of 
the  overall  classification  scheme.  The  process  is  then  repeated  iteratively,  training  a  second 
ML  classifier  using  data  classified  by  the  first.  It  is  shown  that  as  applied  to  polarimetric 
SAR  images  of  San  Francisco  and  Beaufort  sea  ice,  acquired  by  the  Jet  Propulsion  Labo¬ 
ratory,  this  algorithm  converges  rapidly,  and  significantly  improves  classification  accuracy. 
Performance  after  convergence  is  seen  to  be  comparable  to  that  obtained  with  a  supervised 
ML  classifier,  while  maintaining  the  advantages  of  an  unsupervised  technique. 
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For  polarimetric  remote  sensing,  geophysical  m'*^'a  are  modeled  as  layers  containing 
randomly  embedded  scatterers.  In  a  medium  Each  as  vegetation  canopy,  the  scatterers 
can  have  various  shapes,  sizes,  and  perrolUivities  that  are  significantly  distinct  from  the 
background  medium.  The  model  presented  in  this  paper  will  consider  each  type  of  the 
scatterers  as  a  species  which  can  takes  on  a  shape,  size,  and  complex  permittivity  different 
from  other  species.  The  effective  permit*ivity  of  the  random  medium  is  derived  under  the 
strong  permittivity  fluctuation  theory  and  the  polarimetric  scattering  coefficients  are  cal¬ 
culated  for  the  layer  configuration  with  the  distorted  Born  approximation  in  the  analytical 
wave  theory  which  preserves  the  phase  information. 

The  multiple  species  in  the  random  medium  are  considered  as  randomly  oriented  ellip¬ 
soids  and  described  by  multiple  three-dimensional  ellipsoidal  correlation  functions.  The 
variances  and  correlation  lengths  of  the  correlation  functions  characterize  the  fluctuation 
strengths  and  the  physical  geometry  of  the  scatterers,  respectively.  The  strong  permittiv¬ 
ity  fluctuation  theory  is  extended  to  account  for  the  multiple  species  of  ellipsoidal  shape. 
In  the  random  medium,  a  coincidence  of  an  observation  point  with  a  source  point  gives 
rise  to  the  singularity  of  the  dyadic  Green’s  function  which  is  properly  taken  into  account 
with  exclusion  volumes  of  the  scatterers. 

Polarimetric  scattering  properties  of  a  remotely  sensed  medium  are  depicted  with  a 
covariance  matrix  whose  elements  are  polarimetric  scattering  coefficients.  The  medium 
has  a  layer  configuration:  the  top  layer  such  as  air  is  considered  as  homogeneous,  the 
middle  layer  such  as  a  vegetation  canopy  is  random,  and  the  underlying  layer  such  as  soil 
is  a  homogeneous  half  space.  More  random  medium  layers  can  also  be  introduced  in  the 
configuration  to  account  for  weather  effect  such  as  fog  cover.  The  distorted  Born  approx¬ 
imation  is  then  used  with  the  effective  permittivity  to  obtain  the  polarimetric  covariance 
matrix.  The  result  for  the  cross-polarized  return  is  non-zero  even  in  the  first  order 
approximation.  Due  to  the  non-spherical  shape  and  the  random  orientation  of  the  scatter¬ 
ers,  the  correlation  coefficient  between  the  HH  and  VV  returns  has  a  magnitude  different 
from  unity  and  a  small  phase.  The  scattering  coefficients  are  also  used  to  calculate  the 
Mueller  matrix  for  synthesis  of  polarization  signatures.  The  copolarized  signature  of  the 
random  medium  has  a  rather  straight  distortion  track  and  a  recognizable  pedestal. 
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To  demonstrate  the  use  of  polarimetry  in  passive  remote  sensing  of  azimuthal  asymmet¬ 
ric  features  on  a  terrain  surface,  an  experiment  is  designed,  implemented,  and  presented 
in  this  paper. 

The  experiment  setup  is  illustrated  in  the  figure  below.  A  triangular  corrugation  pattern 
with  height  h  =  2.5  cm  and  period  p  —  5  cm  is  made  on  the  sandy  soil  surface  covering 
23  periods  by  length  £  =  160  cm  and  thickness  d  =  12.7  cm.  A  radiometer  of  15°  beam 
width  operating  at  10  GHz  is  used  in  the  experiment.  The  radiometer  is  mounted  on 
a  tripod  at  an  elevation  of  1.8  m  height  and  directed  toward  the  soil  surface  along  the 
look  direction  determined  by  azimuthal  angle  4>  and  polar  angle  9.  Polarimetric  brightness 
temperatures  Tgf^,  and  Tpp  are  meatsured  respectively  with  horizontal,  vertical,  and 

45°  polarization  orientations  for  6  =  20°, 30°  and  (p  from  0°  to  90°.  From  the  measured 
temperatures,  the  third  Stokes  brightness  temperatures  are  obtained. 

Absolute  values  as  high  as  30-40  K  of  the  third  Stokes  brightness  temperatures  are 
observed.  A  theoretical  analysis  of  the  data  indicates  that  the  high  values  oi  Uj]  are 
caused  by  the  azimuthal  asymmetry  on  the  remotely  sensed  soil  surface.  It  is  also  observed 
from  the  experiment  that  decreases,  increases,  T^p  decreases  to  a  minimum  at 

4>  =  45°  and  then  increases  as  <p  takes  on  the  increased  vadues.  For  U the  trend  is  similar 
to  that  of  T'gp.  These  general  trends  are  supported  by  our  theoretical  predictions  of  the 
polarimetric  brightness  temperatures.  The  significant  observation  from  this  experiment 
is  that  the  surface  asymmetry  can  be  detected  with  a  measurement  ol  U ^  at  a  single 
azimuthal  angle. 
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Among  the  various  theoretical  models  applied  to  the  study  of  electromagnetic  wave 
scatterings  from  geophysical  terrain,  such  as  snow,  ice,  and  vegetation  canopy,  the  radiative 
transfer  theory  has  drawn  intensive  attention  in  the  microwave  remote  sensing  society 
during  the  past  years.  In  most  of  the  scattering  models,  the  volume  scattering  and  the 
surface  scattering  effects  have  been  investigated  separately.  Recently,  there  has  been  a 
growing  interest  in  the  construction  of  composite  models  which  can  take  into  account  both 
types  of  scattering. 

In  this  research,  we  present  the  first  order  iterative  solution  to  the  vector  radictive 
transfer  equations  for  a  two-layer  medium  with  a  diffuse  top  boundary  and  an  irregular 
bottom  boundary  of  Gaussian  roughness.  The  Kirchhoff  approximation  and  the  geomet¬ 
rical  optics  approach  with  shadowing  correction  are  used  in  formulating  the  boundary 
conditions.  To  demonstrate  the  utilities  of  this  theory,  we  apply  our  formula  with  a  phase 
matrix  for  randomly  distributed  and  oriented  spheroidad  discrete  scatterers  to  calculate 
the  backscattering  coefficients  from  soybean  field  at  different  growing  stages  and  compare 
the  results  with  experimental  measurements.  Good  agreement  has  been  achieved  for  both 
the  co-polarized  and  the  cross-polarized  data.  It  is  observed  that  the  presence  of  the  rough 
surface  can  significantly  enhance  the  backscattering  at  small  incident  angles  and  the  levels 
of  the  cross-polarized  return.  The  polarization  signatures  calculated  based  on  the  Mueller 
matrix  shows  a  straight  distortion  track  and  an  observable  pedestal.  Numerical  comparison 
to  the  backscattering  coefficients  calculated  by  using  planar  bottom  boundary  conditions 
with  or  without  the  incoherent  addition  of  the  rough  surface  effects  are  also  made  in  this 
work. 
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ABSTRACT.  Polarimetric  calibration  algorithms  using  combinations  of  point  targets 
and  reciprocal  distributed  targets  are  developed.  From  the  reciprocity  relations  of  dis¬ 
tributed  targets,  an  equivalent  point  target  response  is  derived.  Then  the  problem  of  po¬ 
larimetric  caJibration  using  two  point  targets  and  one  distributed  target  reduces  to  that 
using  three  point  targets,  which  has  been  solved  before.  For  calibration  using  one  point 
target  and  one  reciprocal  distributed  target,  two  cases  are  analyzed  with  the  point  tar¬ 
get  being  a  trihedral  reflector  or  a  polarimetric  active  radar  calibrator  (PARC).  For  both 
cases,  the  general  solutions  of  the  system  distortion  matrices  are  written  as  a  product  of 
a  particular  solution  and  a  matrix  with  one  free  parameter.  For  the  trihedral-reflector 
case,  this  free  parameter  is  determined  by  assuming  azimuthal  symmetry  for  the  dis¬ 
tributed  target.  For  the  PARC  case,  knowledge  of  one  ratio  of  two  covariance  matrix 
elements  of  the  distributed  target  is  required  to  solve  for  the  free  parameter.  Numerical 
results  are  simulated  to  demonstrate  the  usefulness  of  the  algorithms  developed  in  this 
paper. 

I.  Introduction 

In  this  paper  the  measured  scattering  matrix  corresponding  to  a  target  with  scat¬ 
tering  matrix  S  is  assumed  to  be  [l] 

X  =  cRST  (1) 

The  matrices  T  and  R  account  for  the  channel  imbalance  and  cross-polarization  coupling 
of  the  transmitting  and  receiving  ports,  respectively.  Propagation  delay  and  loss  between 
the  radar  and  the  target  is  accounted  for  by  parameter  c.  The  objective  of  polarimetric 
calibration  is  to  solve  for  the  R  and  T  matrices  using  responses  from  targets  with  known 
polarization  scattering  parameters. 

For  general  polarimetric  systems,  six  parameters  need  to  be  estimated  for  relative 
polarimetric  calibration;  three  for  each  of  the  transmitting  and  receiving  ports  including 
one  for  channel  imbalance  and  two  for  cross-talk.  To  calibrate  a  general  polarimetric 
system,  algorithms  using  three  point  targets  (in-scenc  reflectors)  have  been  developed  in 
[1-5]  and  a  calibration  algorithm  using  one  distributed  target  and  two  point  targets  was 
presented  in  [6,7|. 

For  polarimetric  radars  with  a  specific  form,  the  number  of  calibration  targets  can 
be  reduced.  For  example,  a  commonly  adopted  model  is  the  reciprocal  radar  {R  being 
the  transpose  of  T).  In  this  case,  a  calibration  technique  using  natural  distributed  tar¬ 
gets  with  azimuthal  symmetry  and  trihedral  corner  reflectors  was  developed  |8]. 

In  this  paper,  we  will  investigate  algorithms  for  calibrating  a  general  polarimetric 
system  using  combinations  of  point  and  distributed  targets. 

II.  Equivalent  Point  Target  Response  Using  Reciprocity 

In  this  section,  we  derive  an  equivalent  point  target  from  the  response  of  reciprocal 
distributed  targets,  utilizing  the  model  of  polarimetric  radars  defined  by  Equation  (1). 
Let  us  define  a  mapping  between  the  scattering  matrix  X  and  a  column  vector  X  as 

=  (2) 
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where  the  superscript  t  signifies  transpose.  Likewise,  we  can  map  the  matrix  5  into  the 
column  vector  S.  Thereafter,  Equation  (1)  can  be  converted  into  a  matrix  equation 

X  =  AS  (3) 

where  A  is  the  direct  product  of  R  and 

A  =  R®T^  (4) 

For  the  case  that  S  is  the  response  from  a  distributed  target,  it  is  more  appropriate 
to  describe  the  distributed  target  by  its  covariance  matrix,  C,.  Prom  Equation  (3),  the 
measured  covariance  matrices  Cx  of  X  is  related  to  Cs  by 

Cx  =  AC,A-^  (5) 

where  *  +  ’  signifies  conjugate  transpose. 

By  assuming  the  distributed  target  is  reciprocal  f5i2  =  "?2l)  covariance 

matrix  Cj  has  rank  3,  Klein  and  Freeman  [6,7]  derivea  the  following  equation 

=  (6) 

where  Qx  and  Q,  are,  respectively,  the  eigenvectors  of  Cx  and  Cm  with  zero  eigenvalue. 

In  addition,  5]2  =  S21  means  that  columns  2  and  3  of  Cm  are  equal.  Hence, 

Qm  =  [Q  -1  1  0)*  (7) 

Note  that  c'  is  an  appropriate  scaling  constant,  taking  into  account  the  arbitrary  abso¬ 
lute  magnitudes  of  eigenvectors. 

Equation  (6)  is  an  important  observation  made  by  Klein  and  Freeman  [6,7]  using 
the  reciprocity  relations  satisfied  by  the  distributed  targets.  Subsequently,  tney  carried 
out  the  matrix  inversion  and  multiplication  in  Equation  (6)  and  obtained  three  nonlinear 
equations  for  matrix  elements  of  R  and  T  [Ref.  6,  Eqs.  (19)-(2l)].  Further  supplement¬ 
ing  these  three  equations  with  the  response  from  two  point  targets,  they  solved  for  the 
normalized  matrix  elements  of  R  and  T  explicitly. 

A.  Equivalent  Point  Target  Response 

In  this  paper,  we  recast  Equation  (6)  into  another  form  to  allow  for  a  general  appli¬ 
cation.  After  realizing  that  the  mapping  from  Equation  (1)  to  Equation  (3)  is  one  to  one 
in  nature,  we  can  convert  Equation  (6)  into  the  following  form 

Qt~^  =  j^RQt~^T  (8) 

where  Q,  and  Qx  are  2  by  2  matrices  corresponding  to  vectors  Qm  and  Gi,  respectively, 
through  the  mapping  defined  by  Equation  (2).  **’  denotes  complex  conjugate.  Compar¬ 
ing  Equation  (8)  with  Equation  (1),  we  can  identify  that  is  the  response  from  a 

point  target  with  scattering  matrix  Qf 

B.  Reciprocation  of  Polarimeiric  Radar 

As  we  know  that  a  polarimetric  radar  is  in  general  not  reciprocal  {R  ^  T^),  there¬ 
fore  the  measured  scattering  matrix  X  is  not  symmetric  when  5  is  symmetric.  As  to  be 
shown  in  the  following.  Equation  (8)  can  be  utilized  to  symmetrize  the  measured  scatter¬ 
ing  matnx. 
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From  Equation  (8),  it  can  be  shown  that 

r  =  i2‘(g+Qi)"^  (9) 

For  convenience,  all  the  proportionaJ  constants  like  c'  will  not  be  written  out  explicitly 
from  now  on.  Then  substituting  Equations  (9)  into  (1),  it  is  easy  to  show  that 

y  =  rsr*  (lo) 

Note  that  if  5  is  symmetric,  so  is  Y.  After  comparing  Equations  (10)  with  (1),  we  can 
readily  interpret  Y  as  the  response  observed  by  a  reciprocal  polarimetric  radar  with  T 
being  the  transpose  of  R. 

If  X  is  the  backscatter  from  a  reciprocal  distributed  target,  then  we  can  derive  from 
C'x  a  reciprocated  covariance  matrix  Cy  of  Y  which  is  given  by 

Cy  =  AyCxAy  (11) 

where 

=  (12) 

The  results  of  this  section  show  that  the  target  reciprocity  can  be  utilized  to  derive 
an  equivalent  target  response  and  makes  a  polarimetric  radar  reciprocal. 

III.  Calibration  Using  Two  Point  Targets  and  One  Reciprocal  Distributed 
Target 

This  section  discusses  how  to  use  two  general  point  targets  and  one  reciprocal  dis¬ 
tributed  target  for  polarimetric  calibration.  Note  that  a  general  class  of  problems  where 
these  two  point  targets  are  two  reciprocal  reflectors  have  been  solved  explicitly  by  Klein 
and  Freeman  and  verified  with  the  Jet  Propulsion  Laboratory  polarimetric  sensors  [7). 
Here,  we  will  show  how  the  equivalent  target  can  be  used  to  solve  the  problem. 

As  shown  in  the  previous  section,  we  obtain  an  equivalent  target  from  the  reci¬ 
procity.  Hence,  by  combining  this  equivalent  target  with  another  two  point  targets,  the 
problem  of  calibration  using  two  point  targets  and  one  distributed  target  is  converted  to 
that  using  three  point  targets,  which  has  ^ready  been  solved  [3). 

IV.  One  Trihedral  Reflector  and  One  Reciprocal  Distributed  Target 

This  section  discusses  how  to  ezdibrate  a  polarimetric  radar  using  the  responses 
from  a  trihedral  reflector  and  one  distributed  target.  Note  that  van  Zyl  [8]  has  already 
provided  a  scheme  for  calibrating  a  reciprocal  polarimetric  radar  (T  =  R})  using  one 
trihedral  reflector  and  one  distributed  target  with  azimuthal  symmetry.  Also  in  the  pre¬ 
vious  section  we  have  shown  that  a  general  polarimetric  radar  (i?  T*)  can  be  made 

symmetric  by  using  the  reciprocity,  indicating  that  van  Zyl’s  approach  can  be  directly 
applied.  However,  we  will  take  a  different  approach  in  solving  the  problem  so  that  par¬ 
tial  calibration  can  be  carried  out  without  assuming  the  azimuthal  symmetry  for  dis¬ 
tributed  targets. 

Suppose  that  the  measured  response  from  a  trihedral  is  given  by 


Xc  =  RT 


(13) 


Hence, 

Kc  =  XcQtQ[  =  RR^ 


(14) 
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and  T  can  be  expressed  as 

T  =  R-^Xc  (15) 

Tlie  general  solution  of  R  can  be  written  as  the  product  of  a  particular  solution  Rp  and 
a  matrix  M 

R  =  RpM-^  (16) 

where  the  particular  solution  Rp  has  the  following  form: 


Rr.= 


^p2\  Rp22 


M  = 


cos  Q  sin  B 
—  sin  B  cos  B 


Note  that  6  is  allowed  to  be  complex.  After  substituting  Equation  (16)  into  Equation 
(14),  we  can  solve  for  the  following  quantities  from  the  resulting  equation 


■^p21  _  l^cl2  _  lc21 

RpW  ^cll  Yell 

^p22  ,  [Yc22 


/V’c22 

ynii 

(19a) 

(195) 


Then  from  Equation  (15),  we  obtain  the  particular  solution  of  Tp  (=  R^^Xc)  and 
the  general  solution  of  T  (=  MTp).  Hence,  the  general  forms  of  R  and  T  have  been 
deduced  and  only  one  parameter  a  (=  tan  6)  remains  to  be  determined.  In  the  follow¬ 
ing,  we  assume  the  azimuthal  symmetry  for  the  distributed  target  to  determine  a.  Az¬ 
imuthal  symmetry  of  the  distributed  target  implies  that 

^^*12  =  CslZ  -  C,24  =  =  0  (20) 

By  substituting  the  general  form  of  R  and  T  into  Equation  (1),  it  can  be  shown  that 

Xp  =  R;^XTp^  =  M-'^SM  (21) 

Note  that  M~^  -  A/‘.  Therefore,  matrix  Xp  is  symmetric. 

Representing  Xp  by  a  vector  Xp  according  to  Equation  (2),  we  can  define  the  co- 
variance  matrix  Cp  =  E{XpX^]  which  is  related  to  Cx  and  C*  by 


Cp  —  XpGxAp 

(22a) 

—  AmO$A^ 

(225) 

Ap  =  Rp^  0  Tf'* 

(23a) 

0  M* 

(235) 

Multiplying  both  sides  of  Equation  (21)  by  Af,  then  calculating  the  ensemble  av¬ 
erage  of  the  products  of  matrix  elements,  and  using  the  fact  that  Cpi2  =  Cpiz  and 
C’p34  =  Cp24i  we  obtain  the  following  equations  for  a 

oiCp22  +  C',22  -  C'p22  -  ^^»32)  +  o*(C'pl4  “  C,ii  +  CpH  -  C,4i)  =  -(C'pl2  +  C'p42)  (24) 
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c^[Cp32  +  <^523  ~  C'p23  -  C'jSa)  +  <3t*(-C'pii  +  Cgi4  —  CpH  +  (7,44)  =  -(^pI3  +  C'p43)  (25) 

Note  that  a  can  be  obtained  by  utilizing  the  solution  of  the  following  equation 

aa  +  ba*  =  c  (26) 

wliich  has  tl^e  following  solution 

Or  =  — » - 5,  provided  la]  ^  \b\.  (27) 

jo|2  _  |i|2’ 

The  criterion  |al  /  j6|  can  be  satisfied  for  a  large  class  of  distributed  targets  as  observed 
from  Equations  (24)  and  (25). 

If  we  start  with  a  polarimetric  radar  having  good  channel  isolation  (small  q),  we 
can  use  the  following  iterative  scheme  to  determine  the  value  of  a: 

Step  1.  Initialize  a  as  zero. 

Step  2.  Evaluate  Am  from  Equation  (23b)  with  the  current  value  of  a  and  caJculate 
Cj  =  A^CpA'^'^  from  Equation  (22b)  as  the  current  estimate  of  C,. 

Step  3.  Solve  a  from  either  Equations  (24)  or  (25)  by  using  the  current  value  of  C'^. 

Step  4.  Repeat  Steps  2  and  3  until  C'jj,  ^*^24'  ^*34  smaller  than  a  se¬ 

lected  threshold,  which  in  our  case  is  set  to  be  10~®  times 
Once  Q  is  solved,  the  solution  for  R  and  T  is  complete  except  for  their  absolute 
magnitudes. 

V.  One  Polarimetric  Active  Radar  Calibrator  and  One  Reciprocal  Distributed 
Target 

This  section  discusses  how  to  use  a  PARC  and  a  distributed  target  for  calibration. 
Suppose  that  there  is  a  PARC  with  the  following  scattering  matrix 

5c  =  2  J  (28) 

Subsequently,  substituting  Equation  (28)  into  (10)  results  in  a  equation  for  the  trans¬ 
formed  scattering  matrix  Yc-  Then  it  is  straightforward  to  show  that  the  general  solution 
of  R  has  the  following  form 

R  =  R^sM-^  (29) 

where  ti-.e  particular  solution  Rp  is 

[ wn., 

and 

Mp  = 

Here  /?  is  a  free  parameter.  In  addition,  the  particular  solution  of  T  is  obtained  from 
Equation  (8)  and  is  given  by 

=  (32) 

It  can  be  shown  that  the  general  solution  of  T  can  be  written  as 

T  =  ^M-%  (33) 


1  0 
0  1//3  • 
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Substituting  Equations  (29)  and  (33)  into  Equation  (1)  results  in 


(34) 


Let  ^  and  vectorize  it  as  by  Equation  (2).  We  can  relate  the  covari¬ 

ance  matrix  =  E[X^X^]  to  Cx  by 


C/i  —  ■A.^iCxA.'^  (35) 

where 

=  R-^  0  (36) 

Multiplying  X^  by  Ry,  from  the  left-hand  side  of  (34)  and  carrying  out  the  ensemble  av¬ 
erage  of  products  of  matrix  elements,  we  can  obtain  the  following  equations 


/  C,ll  Cyl2^ 

• 

\  -  1 

fCsn 

V<7,13 

C'iH  _ 

<7,13 

Cy23 

Cs22  ^7^12 

<7,23 

<7/113 

/ C,22  Oy2\^ 

0 

\  - 

/  <7,23 

<7/134  \ 

\Ct2A  ^/i22y 

\<7,34 

<7/123/ 

-  4.  / ( _  ,  jCsiiCyii 

^  \l\CsHCynJ 


(37a) 

(376) 

(37c) 


(38) 


Here  we  can  see  that  in  order  to  solve  the  parameter  /?,  it  is  necessary  to  know  one  of 
the  ratios  of  matrix  elements  of  C$  beforehand.  Also  notice  that  Equation  (37)  is  not 
applicable  to  distributed  targets  with  azimuthal  symmetry  which  will  result  in  the  case 
of  zero  divided  by  zero. 


VI.  Numerical  Simulation  and  Discussion 

In  this  section  numerical  simulation  of  algorithms  discussed  in  the  previous  two  sec¬ 
tions  will  be  presented.  The  polarization  transfer  matrices  of  the  polarimetric  radar  are 
assumed  to  be 


1  0.0426Z  -  169.5° 

0.0532/113.6°  1.0638/  -  86.3° 


(39) 


1  0.1042^-77.8° 

0.0625^30°  1.0417;! -57.9“ 


(40) 


This  set  of  system  matrices  corresponds  to  the  distortion  matrices  of  the  Jet  Propulsion 
Laboratory  L-band  radar  estimated  by  a  three-point-target  approach  [9).  Note  that  the 
cross-talk  errors  are  smaller  than  -20  dB  and  the  channel  imbalances  are  around  0.5 
dB.  The  most  severe  errors  are  the  differences  in  phases  between  channels,  which  will 
significantly  distort  the  co-polarization  signature  [10]  of  trihedral  reflectors  and  make  it 
look  like  that  of  a  dihedral  reflector. 
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In  this  paper,  we  denote  channels  1  and  2  as  horizontal  (h)  and  vcrticad  (v)  polar¬ 
izations,  respectively.  Therefore,  the  parameters  of  a  covariance  matrix  C$  are  defined 

^  y/^hvPhhhv  yJ^vhPhhvh  \/7Phhvv 

V^hvPhhhv  y/^hv^vhPhvvh  y/^hv'fPhvvv 

y/^vhPhhvh  \/^hv^vhPhvvh  y/^vh'yPvhvv 

y/^Phhvv  '/^hvKPhvvv  y/^vh'7Pvhvv 

Here,  Pa0f6  represents  the  correlation  coefficient  of  ‘a/3’  and  ‘76’  polarizations.  For  dis¬ 
tributed  targets  with  azimuthal  symmetry,  Phhhv  =  Phhvh  =  Phwv  =  Pvhw  =  0- 

Listed  in  column  1,  Table  1  are  the  covariance  parameters  of  the  distributed  target 
to  be  used  for  the  simulation  of  calibration  algorithms,  and  its  co-polarization  signature 
is  shown  in  Figure  1(a).  Due  to  the  channel  imbaJance  and  cross-talk  for  the  radar  de¬ 
fined  by  Equations  f39)  and  (40),  the  observed  co-polarization  signature  would  be  that 
shown  in  Figure  1(b).  After  the  reciprocity  is  applied,  the  co-polarization  signature  cor¬ 
responding  to  the  reciprocated  covariance  matrix  Cy  is  illustrated  in  Figure  1(c),  which 
becomes  symmetric  with  respect  to  the  90-degree  prientation  angle. 

A.  One  Trihedral  Reflector  and  One  Reciprocal  Distributed  Target 

Considered  below  is  the  trihedral-reflector  approach  (see  Section  III).  Two  solutions 
were  obtained.  One  produces  the  correct  R  and  T  and  the  other  is  related  to  the  correct 

one  by 

R'  =  R 

Using  these  two  solutions  to  calibrate  the  radar,  the  calibrated  covariance  parameters  of 
the  distributed  target  are  the  same  as  those  given  in  column  2  of  Table  1.  This  is  due 
to  the  fact  the  sign  difference  does  not  affect  the  responses  from  trihedral  reflectors  and 
azimuthally  symmetric  distributed  targets. 

Particular  solutions  for  the  case  considered  are  given  by 

P  _[  1  0  1  r  -  f  1  0.0698Z  -  96.51°' 

~  [0.098/109.24°  1.06/ -  86.40°  J  ’  “  [0.1038/22.21°  1.0482/ -  57.84° 

and 

p  _[  1  0  1  7,  _  [  1  0.0698/ -  96.51°' 

-  [0.098/109.24°  1.06/93.6°J’  -  |^ 0.1038/ -  157.79°  1.0482/122.16° 

Note  that  the  first  solution  is  very  close  to  the  correct  system  parameters,  in  particular, 
the  channel  imbalance.  If  the  system  is  calibrated  by  the  particular  solution,  the  cali¬ 
brated  distributed  target  parameters  are  given  in  columns  4  and  5  of  Table  1,  and  the 
corresponding  polarization  signatures  are  shown  in  Figures  1(d)  and  1(e),  respectively, 
which  are  not  visibly  distinguishable  from  the  correct  signature  shown  in  Figure  1(a). 

B.  One  Polarimeiric  Active  Radar  Calibrator  and  One  Reciprocal  Distributed  Target 

The  PARC  approach  outlined  in  Section  IV  is  considered  next.  The  particular  solu¬ 
tion  is  given  by 

1  0.400/ -  83.2°  1  r  -  f  1  0.1042/ -  77.8°' 

0.0532/113.6°  1.0  J’  '"“[0.0588/116.3°  0.9792/28.4° 

We  can  see  that  R^2\  ^^d  T^i2  match  correctly  with  R2\  and  T12  [see  (29)  and  (33)|. 
Calibration  by  the  particular  solution  alone  yields  the  parameters  of  0^  which  are  shown 
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in  Tabic  2.  Note  that  all  the  correlation  coefficients  are  well  calibrated  except  the  phase 
Phhvv  Hence,  if  we  have  a  system  with  an  ideal  phzise  balance  between  h  and  v  chan¬ 
nels,  then  all  the  correlation  coefficients  can  be  well  calibrated  by  the  particular  solution, 

and  T;, . 

Further  assuming  that  two  of  the  distributed  target  parameters  7  and  Phhw 
known,  wc  can  apply  Equation  (38)  for  the  calculation  of  /?  and  to  obtain  two  solutions. 
One  correctly  produces  the  system  parameters,  whereas  the  other  is  related  to  it  by 
Equation  (42).  By  calibrating  the  system  using  these  two  solutions,  the  calibrated  co¬ 
variance  matrix  parameters  are  given  in  column  2,  Table  2,  which  correctly  match  the 
original  parameters.  It  should  be  noted  that  if  correlation  coefficients,  Phhhv)  Phhvht 
Phvvv^  s^rid  p^ihvv  zeros,  then  the  correlation  parameters  after  c^ibration  by  the 

wrong  solution  will  carry  an  opposite  sign. 

C.  Comparison  of  Trihedral  Reflector  and  PARC  Approaches 

The  results  of  the  above  two  sub-sections  indicate  that  when  the  particular  solution 
is  applied  for  calibration,  the  trihedral-reflector  approach  is  quite  robust  in  calibrating 
the  channel  imbalance  and  weak  at  detecting  the  cross-polarization  couplings.  In  con¬ 
trast,  the  PARC  approach  is  poor  at  removing  the  channel  imbeJance  and  useful  in  esti¬ 
mating  the  cross-talk.  In  other  words,  the  trihedral-reflector  approach  is  useful  when  the 
system  has  relatively  small  cross-talk,  and  the  PARC  approach  is  useful  when  the  sys¬ 
tem  has  good  channel  balance,  if  no  further  assumption  about  the  distributed  targets  is 
made. 

We  should  ..ote  that  the  alignment  of  point  target  is  less  critical  in  the 
trihedral-reflec' j-  .pproach  than  the  PARC  approach,  since  the  PARCs  needs  to  be 
carefully  orientca  in  order  to  avoid  artificial  cross-talks.  However,  the  PARC  approach 
can  provid  a  larger  signal  to  noise  ratio  than  the  trihedral-reflector  approach. 
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Figure  1.  Co-polarization  signatures  of  a  distributed  target  with  azimuthal  sym¬ 
metry:  (a)  original,  (b)  distorted  (Cx),  (c)  reciprocated  (Cy),  (d)  calibrated  by  the 
particular  solution,  Rp\  and  Tp\,  and  (e)  calibrated  by  the  particular  solution,  Rp2 
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Table  1.  Original,  observed,  and  calibrated  covariance  matrix  parameters  by  using 
the  responses  from  a  trihedral  reflector  and  azimuthal  symmetric  distributed  targets. 
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Table  2.  Original,  observed,  and  calibrated  covariance  matrix  parameters  by  using  the 
responses  from  a  PARC  and  azimuthal  symmetric  distributed  targets  with  known  7  and 
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Abstract-  A  multivariate  A-distribution,  well  supported  by  experimental  data,  is  pro¬ 
posed  to  model  the  statistics  of  fully  polarimetric  radar  clutter  of  earth  terrain.  In  this 
approach,  correlated  polarizations  of  backscattered  radar  returns  are  characterized  by  a 
covariance  matrix  and  homogeneity  of  terrain  scatterers  is  characterized  by  a  parameter 
a.  As  compared  with  C-,  L-  and  P-band  polarimetric  SAR  image  simultaneously  mea¬ 
sured  by  Jet  Propulsion  Laboratory  (JPL)  on  Mt.  Shasta,  it  is  found  that  a  appears  to 
decrease  &om  C-  to  P-band  for  both  the  forest  and  burned  areas. 

I.  INTRODUCTION 

Polarimetric  terrain  backscatter  data  observed  with  satellite  and  airborne  syn¬ 
thetic  aperture  radars  (SAR)  have  demonstrated  potential  applications  in  geologic 
mapping  and  terrain  cover  classification  [1-7|.  In  previous  publications  on  this 
subject,  Gaussian  statistics  have  been  frequently  assumed  for  the  radu  return 
signals'to  build  the  Bayes  terrain  classifier  [1-3].  However,  abundant  experimen¬ 
tal  evidence  shows  that  terrain  radar  clutter  is  non-Gaussian,  i.e.,  non-Rayleigh 
in  amplitude  distribution  [8-10]. 

Among  many  non-Gaussian  statistics,  the  if -distribution  has  proven  to  be 
j  useful  in  characterizing  the  amplitude  distribution  of  electromagnetic  echoes  from 
various  objects  [9-15],  including  diverse  ground  surfaces  [9],  sea  surface  [10]  and 
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CALIBRATION  OF  POLARIMETRIC  RADARS 
USING  IN-SCENE  REFLECTORS 


S.  H.  Yueh,  J.  A.  Kong,  and  R.  T.  Skin 


9.1  Introduction 

9.2  Polarimetric  Calibration  Using  Three  In-Scene 
Reflectors 

a.  Si  and  5:  are  Singular 

b.  Si  and  S2  are  Nonsingulaur 

c.  Si  is  Singular  and  S2  is  Nonsingular 

9.3  Uniqueness 

9.4  Polarimetric  Calibration  Using  Three  Simple  In-Scene 
Reflectors 

a.  Two  Dipoles  and  One  General  Target 

b.  One  Trihedral  Reflector,  One  Target  with  Diagonal  Scattering 
Matrix,  and  One  General  Target 

c.  One  Dipole,  One  'li’ihedral  Reflector,  and  One  General  Reflector 

9.5  Misalignment 

a.  Uniform  Misalignment 

b.  Relative  Misaligiunent 

9.6  Eflfect  of  Noise 

9.7  Polarimetric  Compensation 

9.8  Summary 
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shown  that  there  are  only  two  linear  independent  roll-inv£iri{knt  targets. 
All  other  roll-invariant  targets  must  be  a  linear  combination  of  these 
two  roll-invariant  targets. 
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CLASSIFICATION  OF  EARTH  TERRAIN 
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6.1  Introduction 

Classification  of  earth  terrun  within  an  image  is  one  of  the  many 
important  applications  of  polarimetric  data.  A  systematic  classification 
procedure  will  place  the  classification  process  on  a  more  quantitative 
level  and  reduce  the  amount  of  photo-interpretation  necessary  (Frons, 
1986].  Single  feature  and  multifrequency  classifications  have  been  used 
in  the  past,  but  classification  can  now  be  applied  to  fully  polarimet¬ 
ric  data  which  have  become  available  due  to  recent  developments  in 
radar  technology.  It  has  been  shown  that  Bayes  classification  using 
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Instead  of  using  the  fully  polarimetric  information,  specific  fea¬ 
tures  can  be  isolated  for  classification  by  deriving  the  appropriate  prob¬ 
ability  density  functions.  The  results  will  then  indicate  the  strong  and 
weak  points  of  each  feature  in  the  classification  of  earth  terredn.  To 
derive  the  classification  procedures  for  simple  features,  the  following 
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Scattering  from  Randomly  Oriented  Scatterers  with  Strong 
Permittivity  Fluctuations 
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R.  T.  Shin 
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Lexington,  MA  02173,  USA 

'  Abstract-  Strong  permittivity  fluctuation  theory  is  used  to  solve  the  problem  of  scatter¬ 
ing  from  a  medium  composed  of  completely  randomly  oriented  scatterers  under  the  low 
frequency  limit.  Based  on  Finkel'berg’s  approach  [2,3],  Gaussian  statistics  is  not  assumed 
for  the  renormalized  scattering  sources.  The  effective  permittivity  is  obtained  under  the 
low  frequency  limit  and  the  result  is  shown  to  be  isotropic  due  to  no  preferred  direction 

i  in  the  orientation  of  the  scatterers.  Numerical  results  of  the  effective  permittivity  are 
illustrated  for  oblate  and  prolate  spheroidal  scatterers  and  compared  with  the  results  for 
spherical  scatterers.  The  results  derived  are  shown  to  be  consistent  with  the  discrete  scat- 
terer  theory.  The  effective  permittivity  of  random  medium  embedded  with  nonspherical 
scatterers  shows  a  higher  imaginary  part  than  that  of  spherical  scatterer  case  with  equal 
correlation  volume.  Under  the  distorted  Born  approximation,  the  polarimetric  covariance 
matrix  for  the  backscattered  electric  field  is  c^culated  for  the  half-space  randomly  ori¬ 
ented  scatterers.  The  nonspherical  geometry  of  the  scatterers  shows  signifleemt  effects  on 
the  cross-polarized  backscattering  returns  os,  and  the  correlation  coefficient  p  between 
HH  and  VV  returns.  The  polarimetric  backscattering  scattering  coefficients  can  provide 
useful  information  in  distinguishing  the  geometry  of  scatterers. 


I.  INTRODUCTION 

There  has  been  a  continually  growing  interest  in  scattering  from  random  media 
embedded  with  nonspherical  particles  because  the  particles  in  rain,  ice  crystals, 
fog,  snow,  leaves,  etc.,  arc  nonspherical.  The  radiative  transfer  theory  has  been 
used  to  study  a  layer  of  randomly  positioned  and  oriented  nonspherical  particles 
overlying  a  homogeneous  dielectric  half  space  [17,15].  A  second  approach  based  on 
wave  theory  and  Foldy’s  approximation  [6]  ,  which  is  valid  only  for  sparse  concen¬ 
trations  of  scatterers  [5]  has  also  been  developed.  For  nonsparse  concentrations  of 
scatterers,  incorporating  the  correlation  of  scatterer  positions,  the  quasicrystalline 
approximation  and  distorted  Born  approximation  have  been  used  to  calculate  the 
mean  field  and  the  incoherent  field  [12,14]. 

I  For  a  medium  with  large  spatial  permittivity  variations,  such  as  water  mois- 

I 

» 
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ture  in  vegetation  and  brine  inclusions  sea  ice,  strong  fluctuation  theory  has  been 
developed  for  uniformly  aligned  spheroidal  scatterers  [11,13,16].  Because  the  sin¬ 
gularity  of  the  dyadic  Green’s  function  is  properly  taken  into  account  in  the  renor¬ 
malization  method,  the  theory  is  applicable  to  both  small  and  large  variances  of 
permittivity  function.  However,  the  orientations  of  scatterers  for  natural  terrain 
media  are  in  general  random,  and,  in  this  paper,  we  extend  the  theory  to  take 
into  account  the  effect  of  random  orientations. 

In  the  derivation  of  strong  fluctuation  theory,  a  Gaussian  random  process  is 
usually  assumed  for  the  renormalized  scattering  sources  [9,10,11,16].  The  Gaus¬ 
sian  statistics  assumption  is  employed  to  factor  the  higher  order  moments  of  the 
fluctuations  in  terms  of  the  products  of  second  moments,  and  a  Feynman  diagram¬ 
matic  representation  of  Dyson’s  equation  for  the  mean  field  is  derived.  Called  the 
bilocal  approximation,  the  mass  operator  in  the  form  of  an  infinite  series  is  ap¬ 
proximated  by  the  first  term.  The  effective  permittivity  is  derived  through  the 
constitutive  relation  between  the  mean  electric  field  and  the  mean  dielectric  dis¬ 
placement.  In  the  case  of  non-Gaussian  fluctuations,  Finkel’berg  [2,3,4]  introduced 
the  one-group  approximation  for  the  mass  operator.  When  the  fluctuations  are 
Gaussian,  the  one-group  approximation  reduces  to  the  Gaussian  case.  In  this 
paper,  Finkel’berg’s  approach  is  generalized  to  randomly  oriented  nonspherical 
scatterers.  The  polarizability  operator  is  also  in  the  form  of  an  infinite  series. 
The  first  term  of  the  series  is  the  same  for  either  Gaussian  or  non-Gaussian  fluc¬ 
tuations.  This  indicates  that  if  the  series  is  approximated  by  the  first  term,  the 
bilocal  approximation  for  the  effective  permittivity  gives  rise  to  the  same  result 
for  Gaussian  and  non-Gaussian  fluctuations. 

In  Section  II,  the  effective  permittivity  for  an  unbounded  medium  embedded 
with  isotropically  oriented  nonspherical  scatterers  is  derived  using  the  strong  fluc¬ 
tuation  theory.  The  effective  permittivity  tensor  is  expressed  in  terms  of  an  infinite 
series.  The  mean  field  is  shown  to  propagate  inside  a  medium  with  the  effective 
permittivity.  The  series  is  truncated  to  the  lowest-order  correction  term  under 
the  bilocal  approximation.  An  analytical  expression  for  the  effective  permittivity 
is  given  for  a  specific  correlation  function  under  the  low  frequency  approximation. 
In  Section  III,  numerical  results  for  the  effective  permittivities  for  spheroidal  scat¬ 
terers  are  compared  with  the  discrete  scatterer  theory  in  the  common  regime  of 
validity.  In  Section  IV,  under  the  distorted  Born  approximation,  the  scattering 
from  a  half-space  with  ellipsoidal  scatterers  is  solved.  The  covariance  matrix  of 
the  polarimetric  backscattering  coefficients  is  obtained.  Because  of  the  azimuthal 
symmetry  of  the  orientation,  four  out  of  nine  elements  of  the  covariance  matrix 
are  shown  to  be  zero  [1].  In  Section  V,  the  polarimetric  backscattering  coefficients 
are  illustrated  for  spheroidal  scatterers  with  different  sixial  lengths. 

II.  STRONG  FLUCTUATION  THEORY 

Consider  an  unbounded  space  filled  with  randomly  distributed  and  oriented  scat¬ 
terers.  Assuming  the  medium  is  characterized  by  the  permittivity  c(r),  the  wave 
equation  for  a  time  harmonic  electric  field  £(r)  is 


I 
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V  X  V  X  E{f)  -kV-^  E(r)  =  0 
«0 


(1) 


where 


f  <5,  r  E  scatterer 
1  ffc,  r  E  background 


(2) 


Introducing  tj  which  will  turn  out  to  be  the  lowest  order  term  of  the  effective 
permittivity  by  solving  (18),  we  can  recast  (1)  into 


V  X  V  X  E(r)  -  fcg  ^  E(r)  =  kl  ^  E{r)  (3) 

«0  <0 

By  introducing  the  dyadic  Green’s  function 

V  X  V  X  dg(fy)  -  kl  ^  dgir,^)  =  Js{r  -  r')  (4) 

(3)  is  transformed  into  the  following  integral  equation: 

E(r)  =Ea  +  klJdT‘  dgir^) .  ~  -)  •  W)  (5) 

Eg  is  the  applied  electric  field.  In  consideration  of  the  singularity  at  the  source 
point,  the  dyadic  Green’s  function  is  decomposed  into  two  parts: 

5g(f,r')  =  PS^5(r,f')-^6(r-r')  (6) 

*0 

where  PS  denotes  the  principal  value  and  S(r)  takes  into  account  the  shape  of 
exclusion  volume  at  r  [16]. 

By  substituting  (6)  into  (5),  the  following  equation  is  obtained: 


Fif)  =  Eo  +  kl  I  df'PSGgiT,r‘)-ar)-Fir') 


(7) 


where  F{r)  is  the  external  field  and  <o((^)  ’  F{f)  corresponds  to  the  induced 
polarization  density: 


F{r)  = 

l{r)  = 


(8a) 

(86) 


Expand  the  solution  of  the  integral  equation  (7)  in  terms  of  infinite  series, 


F(f)  =  £o  +  *.'0  J Gj(r,ri)  •  a?l)  ■  -Eo(n) 

+  ^0  y  y  (P'zFS  Gg{f,ri)  ■  f(ri)  •  PS  Gg(ri,r2)  •  f(r2)  '  •£^0(^2) 


. . 
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+  k^  Jdrijdf2  Jdr^PS  Gg{r,ri}-({fi)  ■  PS  Gg{ri,r2)  ■  ((r2) 

■  PSGg(f2,r3)  -((rj)  •  £0(^3) 

+  •••  (9) 

The  above  equation  can  be  represented  in  an  operator  form 
F(r)  =  Eq  +PS  Gg(r, ri)  ■  ^(rj)  Fo(^l) 

+  PSGg(r,ri).((r:)-PSGgin,r2)-i(T2) 

+  FF(7j(r,fi)  •  ^(rj)  •  PS  Gg(ri,r2)  ■  f(r2)  •  PS  Gg(r2,r3)  ■  ((r^)  ■  Eoirj) 

+  •••  (10) 

Taking  the  ensemble  average  of  (10),  we  obtain 

(F)  =  (7  +  (PSGg  ■T)  +  {PSGg-i.  PSdg  ■  j) 

+  {PSdg-(-PSdg-j  PSdg  T^  +■■■')  -Eq  (11) 

Multiplying  (10)  by  ((r)  and  taking  the  ensemble  average, 

(!•  ^)  =  ((!)  +  (!  •  PS'&g  •!)  +  (!•  PS^g  ■(■PSGg-i) 

+  Q-PS%-tPSVg-hPSGg-j)+--)  -Eo  (12) 

Through  cancelling  Eq  from  (12)  by  using  (11),  the  following  equation  is  obtained; 

(!•?)  =  to!*// (F>  (13) 

where 

Iff  =  (!)  +  ((!  •  PS^,-j)  -  (!)  •  {PSdg.j)) 

+  {{tPSdg-1-PSGg-i)-{tPSdg-t)-{PSdg  l) 

-{j)-{PS^g-i.PSGg-j)+{j)-{PSGg-j)-{PSdg-i)) 

+  •••  (14) 

Note  that  (13)  relates  the  mean  induced  polarization  density  to  the  mean  exter¬ 
nal  field.  Making  use  of  (8),  we  can  relate  the  mean  induced  polarization  density 
and  the  mean  external  field  to  the  mean  dielectric  displacement  (D)  and  the  mean 
electric  field  (£)  .  Thus  the  following  constitutive  relation  can  be  obtained: 

(D{r))  =  <F(r))  (15) 

h/f  =  -f  to  (7  -  !*//  (f ))  !<// 


where 


(16) 
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is  the  effective  permittivity  tensor  of  the  random  medium.  Note  that  in_the 
derivation  of  (15)  and  (16),  the  fact  that  the  D  and  E  are  independent  of  5  or 
the  shape  of  exclusion  volume  has  been  used  [18]. 

By  taking  the  ensemble  average  of  (3)  and  making  use  of  (15),  the  mean  field 
can  be  shown  to  satisfy  the  following  wave  equation: 

V  X  V  X  (£)  -  Jfcg  (E)  =  0  (17) 

*0 

This  indicates  that  the  mean  field  inside  the  random  medium  propagates  in  a 
homogeneous  medium  with  the  effective  permittivity  tensor  ?«//• 


Effective  Permittivity  Under  Low  Frequency  Approximation 

In  evaluating  the  effective  permittivity  (16),  (g  and  fe//  J^ced  to  be  computed. 
Note  that  gives  the  lowest  order  correction  to  the  effective  permittivity  as 

indicated  by  (l4)  and  (16).  In  this  paper,  tg  will  be  adjusted  so  that  the  first 
term  of  igfj,  (14),  is  zero,  i.e., 

(!)  =  0  (18) 

Also,  the  infinite  series  (14)  will  be  approximated  by  the  lowest-order  term 

tffC^{(  PSWg-^)  (19) 

The  above  approximations  correspond  to  the  bilocal  approximation  for  strong 
fluctuation  theory  [17j. 

Assuming 


<  1 


(20) 


the  effective  permittivity  of  the  medium  under  the  low  frequency  approximation 
is  given  by  [Appendix  A] 

=  =(0) 

«e//  -  + «0^e//  (21) 


=  (0) 


I 


(22) 


where  is  given  in  (A6).  Explicitly, 

t”’/  -  jded4>P{e[Ti<t>{T)) 

J  -  +  I(r)  •?(?)  •f(r)j^(r),<^(f)^ 

i  where  0{f)  and  4>{t)  represent  the  corresponding  orientation  angle  of  the  scat- 
i  lerer  at  f  and  P{8,  4>)  corresponds  to  the  probability  density  function  of  the 
!  orientation, 
j  Let 


4(^) 


^(") 


m; 


(23) 
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and  assume  the  normalized  correlation  function  for  the  cross-correlation  between 
the  components  of  ((r)  and  ^(t^)  to  be 


r>  ! - -  ilmj\B,4>  V 


Using  Einstein’s  convention 


(kl  J  |(r)  •  dg{r,^)  ■  iif')  dr"  f(r)  •  f  (F)  •  !{F)|^(r ),  ,^(F)^ 

"  (25) 

=  “^0) 

S  will  be  chosen  in  such  a  way  that  approaches  zero  when  the  frequency 

goes  to  zero. 

In  this  manner  tg  corresponds  to  the  effective  permittivity  at  very  low  frequency. 
Assuming  the  scatterers  are  of  spheroidal  shape,  let 

.  f  r  —  1  .r.fz>ooi 


A:oy  iZ9=o.,^=.0  =  Q  ® 


In  0  0  ] 
0  7p  0  I 


=  r  5p  0  0  I 

5/  =  0  Sp  0  =  —  lim 


U 

0  5, 


1  U  Jo  U  I 

no  6  7,j 


Defining  the  local  coordinates  with  the  following  transformation  [Fig.  1] 
63  =  z  sin  6  cos  ^  -1-  y  sin  S  sin  -f  z  cos  $ 

ej  =  i  sin  ^  —  y  cos  4> 

€2  =  i  cos  9  cos  d’  -f  y  cos  9  sin  4>  —  i  sin  9 
such  that  63  is  in  the  direction  of  the  ajcis  of  the  scatterer.  Therefore, 

j  =  f~^  -I -f 
f  = 

where  the  effective  scattering  source  in  the  scatterer  coordinates  is 
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and  the  coordinate  transformation  matrix  is  given  by 


Figure  1.  Coordinate  system  and  scatterer  geometry. 

Because  the  scatterers  are  assumed  to  be  isotropically  oriented,  the  probability 
density  function  of  the  orientation  is 

P(0,0)  =  ^  (33) 

After  straightforward  algebraic  manipulation,  the  following  result  is  obtained: 

-  [/  (5  U  +  \  f.z)  +(!-/)  Q  ibp  +  \  (fc.)]  7  (34) 

and 
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where 


=  {/  [I  d  (^P  +  ^P)  +  1  «« 
+(1-/) 


(•P  - 


*0  +  ^pi^f  ~  ^g) 


(35) 

{36a) 


= 

ffcp  ~ 


ts  ~  f-g 


*0  +  ^z{*J  ~  ^ff) 
-  tg 

to  +  ^pi^b  ~  *9) 


(366) 

(36c) 


ibz  = 


tfc  -  tg 

to  +  ■5z(ti,  —  tg) 


{36d) 


and  /  is  the  fractional  volume  of  the  scatterers.  The  subscripts  s  and  6  represent 
scatterer  and  background,  respectively. 

For  a  given  correlation  function  (24),  e.g.,  (37),  Sp  and  Sx  are  related  to  €g 
by  (28).  The  equation  for  Cg  is  obtained  by  letting  (34)  equal  to  zero  [see  (18)] 
for  given  Sp  and  Sz-  These  three  coupled  equations  constitute  the  basis  for  the 
solution  of  eg,  Sp,  and  Sz- 

In  this  paper  the  correlation  function  in  the  local  coordinates  (61,62,63)  is 
chosen  to  be 

=  +  (37) 


For  this  correlation  function,  analytical  expressions  of  Sp,  Sz,  Jp,  and  Iz  are 
given  in  Appendix  B. 


m.  NUMERICAL  RESULTS  OF  EFFECTIVE  PERMITTIVITY 


In  this  section,  we  compare  the  effective  permittivities  obtained  using  strong  fluc¬ 
tuation  theory  (21)  and  discrete  scatterer  theory  (Cl)  for  randomly  oriented 
spheroids  under  the  low  frequency  and  small  fractional  volume  limits.  The  effec¬ 
tive  permittivity  obtained  with  the  discrete  scatterer  theory  is  given  in  Appendix 
C.  The  effective  permittivities  are  plotted  as  a  function  of  frequency.  The  parame¬ 
ters  of  the  correlation  function  (37)  are  related  to  the  discrete  scatterer  parameters 
in  the  following  manner.  The  correlation  volume  is  selected  to  be  the  volume  of 
the  discrete  scatterers,  i.e., 


47r  2  .  ^TT  2 


(38) 


and  the  equi-correlation  surface  is  chosen  according  to  the  shape  of  the  scatterer 


\p. 

Iz 


a 

c 


(39) 


Imaflnary  pari  of  e 
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Fi'^urc  2  Imagiuciry  pari  of  as  a  function  of  frequency,  (a)  Ip  = 

•1  min  ai.;l  /;  =  0.0G2.r>inm,  (b)  Ip  =  2mm  ami  1;  =  0.26  mm, 
(c)  Ip  =■  1  inm  ami  Iz  —  1mm,  (d)  Ip  —  0.6mm  and  /;  = 
4  mm,  and  (c)  Ip  =  0.26  mm  and  /.  =  16  mm.  Tlie  daslied 
line  corrcsijonds  to  discrete  scattcrer  llieory.  The  solid  curve 
corresponds  to  strong  fluctuation  theory. 


R*«l  p.rt  «f  lma|ln»ry  part  or  t 
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(a)  l,  •  4  mm,  *  0.0025  mm 


(b)  If  *  2  mm.  I«  •  0.25  mm 
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Figure  3.  Real  part  of  as  a  function  of  frequency,  (a)  Ip  —  4  mm  and 
Iz  =  0.0625  mm,  (b)  Ip  =  2  mm  and  Iz  =  0.25  mm,  (c)  Ip  — 
1  mm  and  I*  =  1  mm,  (d)  Ip  =  0.5  mm  and  Iz  =  4  mm,  and  (c) 
Ip  =  0.25  mm  and  Iz  =  16  mm.  The  dashed  line  corresponds  to 
discrete  scattercr  theory.  The  solid  curve  corresponds  to  strong 
fluctuation  theory. 
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rrw|u«nc7  (OH.)  Fr<qu<nc7  (OHi) 

Figure  3.  Continued. 

The  imaginary  and  real  parts  of  the  effective  permittivities  arc  shown  in  Figs.  2 
and  3,  respectively.  The  permittivities  of  the  scatterers  and  background  are  as¬ 
sumed  to  be  real  so  that  the  scattering  will  be  the  dominant  effect  compared  to  the 
absorption  effect.  The  fractional  volume  of  the  scatterers  is  0.5%.  The  correlation 
volume  (38)  is  kept  a  constant  for  all  cases.  Figures  2(a)  and  2(b)  correspond  to 
disk-like  scatterers  with  correlation  lengths  {Ip  =  4.0  mm  and  /j  =  0.0625  mm) 
and  [Ip  =  2.0  mm  and  =  0.25  mm),  respectively.  Figures  2(d)  and  2(e)  are  for 
needle-like  scatterers  with  correlation  lengths  {Ip  =  0.5  mm  and  Iz  =  4.0  mm), 
and  {Ip  =  0.25  mm  and  Iz  =  16  mm),  respectively.  Figure  2(c)  is  for  spherical 
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scatterers  (/p  =  1.0  mm  and  =  1.0  mm).  The  corresponding  cases  for  the 
real  parts  of  are  shown  in  Figs.  3(a)  through  3(e).  Note  that  the  real  and 

imaginary  parts  of  for  the  nonspherical  scatterers  are  larger  than  the  spher¬ 
ical  scatterer  case.  This  indicates  that  the  nonspherical  scatterer  has  a  higher 
scattering  effect  than  the  spherical  scatterer.  Also,  the  effective  permittivities 
obtained  using  strong  fluctuation  theory  and  discrete  scatterer  theory  are  in  good 
agreement  with  each  other.  The  noticeable  discrepancies  in  Fig.  2(e)  for  frequen¬ 
cies  greater  than  3  GHz  (which  makes  klz  >  1)  is  due  to  the  fact  that  the  low 
frequency  approximation  is  used  to  derive  results  for  both  approaches. 

IV.  BACKSCATTERING  FROM  A  HALF  SPACE  RANDOM  MEDIUM  ^ 

Considering  an  incident  electric  field  impinging  on  a  half  space  filled  with  ran¬ 
domly  distributed  and  oriented  ellipsoid  scatterers  [Fig.  4): 

Eoi  =  jA(-A:ozi)i^Ai  +  v(-l:oz,)£^vt]  «p  {ikpi  ■  p  -  ikjiz)  (40) 

The  unperturbed  field  inside  the  medium  is 

=  [TA(^i.)fil(-fclii)-E^A,-  +  T»(l:i,)ri(-fcui)-E^vi]  exp  (ik^^  ■  p  -  ik^^^z) 

(41) 

where 

2k^ 


Th{k^)  = 
Tv(fcz)= 


kz  -f  fcjx 
2<i  fcz 


^0 


The  scattered  field  in  the  far-field  is  given  as 

foo  rO 


=  kl  r  dp,  [  driGoi(r,r,)  •  ^(r,)  •  F(r,) 
J^OO  •'-oo 


(42a) 

(42i) 

(43) 


where  Gqi  is  the  dyadic  Green’s  function  given  by  [16] 


Goi(?y)  = 


gikor 

47rr 


Xoi{kz)h{kz)h,{ku)+^Yo,{kz)v{kz)v,{ku) 


where 


exp{-ik,  -rj) 

(44) 

A’oi  = 

2kz 

(45a) 

kz  +  kit 

Foi  = 

2tikz 

(456) 

t\kz  + 

h(kz)= 

iky  -  ykx 

kp 

(46a) 

v{kz)  = 

—  zkxkz  —  i/kykz  +  zkp 

(466) 

kokp 
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L  lU  \ 

kp 

(46c) 

.  ,  -xkxku  -  ykyki^  +  ikj 

= - r-7 - - 

kikp 

(46<i) 

kx  =  kfj  cos  4>,  sin  $s 

(47a) 

ky  =  A:o  sin  4>,  sin  6^ 

(476) 

kz  =  ko  cos  dz 

(47c) 

kxi  =  ^0  cos^i  sin  6i 

(47d) 

kyi  =  fcosin^i  sintfj 

(47c) 

^zi  =  ^0  cos 

(47/) 

Z 


Figure  4.  Configuration  for  scattering  from  haJf-space  with  randomly  ori¬ 
ented  spheroids. 

Applying  the  distorted  Born  approximation  to  calculating  the  backscattering 
coefRcients,  we  let 


F(r) 


in 


(48) 


Therefore,  by  substituting  (48)  into  (43),  the  scattered  field  is  given  as 

2  coo  cO  _  _ 

|A'oifi(^-i)  fii(kii)  +  ^  Tbl^'(^i)  vi(A:ij)j 

•  f(^l)  •  [^/i^l(-*l2i)-£^Ai  +  T’wt’l(-fclii)^vi] 


(49) 
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The  scattered  field  is  then  expressed  in  the  following  form: 

’  ^hs  _  fhh  fhv  ^hi 
.  Ev!  ,  ^  .  /vH  fvv  .  .  Eyi 

where  E^^^  and  Evs  are  the  horizontally  and  vertically  polarized  electric  'lelds 
respectively. 

Note  that  for  backscattering  from  a  reciprocad  medium  equals  to  Th' 
polarimetric  covariance  matrix  in  the  backscattering  direction  is  denoted  as  [1] 


<^hh  t^hhhv  <^hh' 

yt)  1 

^hhhv  ‘^kvvv 

^hfivv  ^hvvv 

If 

ifhHfL) 

(Ah/«*v) 

(lAvl") 

{fhvf'.v) 

ifhvn.) 

(l/vvl^) 

Let  the  spectral  density  function  be 

j  ^  cxp(ik  ■  r ) 

The  backscattering  scattering  coefficients  are  given  by 

\Xoin\^  J^’d4>  £de  ^ 

|/j6p  +  -  (tp)  sin^  ^  sin^(9ij  -  (^)j 

+  (1  -/)|efrp  +  Uir-ap)  sin^O  sin2(,^i-,A)f 

=rrl:^iyoirvlH^f  r d4>  C dO^^ 

I  1  ^0  JO  47r 


/  Gp  +  (Or  -  itp)  (  sin  0  cos{<pi  -  d>)+  coiff 
\  «1  *1 


X 

2  ^ 

+  (!-/)  ibp  +  ((bi  -  (bp)  sin  0  cos{4>i  -  <t>)  +  -^  cos 


9  {2kpi,e,<i>) 


..  r  ..sintf 


<rhv  =^4\XoiT„\^  [  rd0 

Jo  Jo 


'  "-^cr  v;ar»«< 
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I  +  k6x-^6p|^] 


sin^  0  sin^(i^,  —  4>) 


sin^  cos(^i  —  ^)  +  COS0 
ki  fc] 


J 


^  {2kpi,e,<j>) 

(53c) 


and  the  cross  correlation  between  polarizations 

•  sine 


Xo  1 7X  >0*1  ^v*  ( r )  I 

\kij  Jo  Jo  4ir 


I  /  [^Jp  +  iisz  -  (sfi)  sin^  8  sin^(<^i  -  <^)] 

^  6p  +  -  6p)  ^  cos((pi  "  ^)  +  "j^ 

+  (>-/)  [ibp  +  Ubz  -  (bp)  ^  sia^i^i  ~  (^)| 

(bp  +  i(bz  -  (bp)  sin0cos{4i  -  ^)  +  ^  cose^ 

9  (2k pi,  0,^) 


2^ 


('hhhv  =  0 
‘^kvvv  ~  ® 


(54a) 

(545) 

(54c) 


where 


-  /■**  1  _ 

<!/{2kp„8,,i>)=  dkz  +  2kui){kz  +  2kl,) 

Under  the  low  frequency  approximation, 

<i{2kp„0,d)  =  ^<^g,p(O)  (56) 

2^1  r, 

Thus 

<^Hb  =  ^iA'oiTJ^  I  lup  +  ^|(f„  -  ^.p)P  +  U;p(6x  -  (,p)]  I 
+  (!-/){  k6pi'  +  l\i(bz  -  ^6p)|"  +  ^Rc  [4‘p(^6x  -  ^6p)]  }  ]  ^e,^{0) 


(57a) 


TT^fen  , 

=  1777^  l>0ir.|* 


2k 


III 


8  l^-p.|^<7”L. 

15 


f  \  |^»p|  +  j  ((^*X  —  (»p)l  +  jR«  \('$p{(zz  —  {»p)] 

|6x  -  ^,p|'|  +(!-/){  \(bp?  +  k^frz  -  ^6p)|' 
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+  jR*  -  ffcp)]  +  ^  -  ^»pl'})  ^^.^(0)  (576) 

f\{($z  -  ^<p)| 


lY  T  1^  ^  l^l*«l  I  ,1/^  ,  \|2 

'■"  5?“'  “■  •'  15 — 


"lit 


+  (l-/)|(az-e6p)|']^8,<t(0) 

<^hhvv  =  ^  XoiT^Yo\T:  (^)  (f\ 

2fcj^.  \Kl/  V  [ 

+  iu«,. -£.,)•  [(^)%(^)'' 


I^jp|  +  T^Spi^lZ  ~  ^>p) 


+ 


+  (1 


+ 


-  /)|  IfipP  +  ^dp  {(bz  -  (bp)  +  ^^6p  ((bi  -  (bp)* 

i|£..-«U’[(^)%(^)]'))»«(0) 


(57c) 


*1 


(57d) 


V.  NUMERICAL  EXAMPLES  OF  POLARIMETRJC  COVARIANCE  MATRICES 


In  this  section  the  results  of  Section  IV  are  illustrated  using  the  parameters  ap¬ 
plicable  to  microwave  remote  sensing  of  vegetation  and  forests.  The  orientations 
of  nonspherical  scatterers  are  assumed  to  be  completely  random.  The  correlation 
function  (37)  is  used  to  evaluate  all  of  the  numerical  results.  Because  four  elements 
of  the  covariance  matrix  are  zero,  it  is  more  convenient  to  write  the  polarimetric 
covariance  matrix  (51)  into  the  following  form  [1] 


C  =  <^hh 


1 

0 

P*^/7 


0 

e 

0 


Pn/t" 

0 

7 


(58) 


where  <  and  7  correspond  to  the  ratios  and  crvvl<^h.h<  respectively,  and 

p  is  the  complex  correlation  coefficient  of  of  HH  and  VV  polarization  returns. 
The  polarimetric  covariance  matrix  contains  very  important  information  in  the 
classification  of  various  terrain  elements  [7]. 

The  effects  of  nonspherical  geometry  on  polarimetric  covariance  matrices  arc 
shown  in  Tables  1  and  2.  Included  are  results  for  five  types  of  scatterers  discussed 
in  Section  III.  Note  that  af^)^  for  the  sphericzJ  scatterer  case  is  bigger  than  the 
nonspherical  scatterer  cues.  This  is  due  to  the  fact  that  the  scattering  loss  (imag¬ 
inary  part  of  <(//)  fof  the  spherical  scatterer  case  is  the  smallest  among  all  cases. 
The  wave  can  penetrate  deeper  into  the  random  medium  with  smaller  scattering 
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loss.  As  a  result,  for  a  haJf-space  configuration,  the  total  scattering  increases.  Also 
note  that  the  ratio  crvvl<^hh  close  to  1  for  all  cases.  This  is  because  all 

the  scatterers  are  completely  randomly  oriented  so  that  the  number  of  scatterers 
which  contribute  to  the  HH  return  is  the  same  as  the  number  of  scatterers  for 
the  VV  return.  The  reason  why  the  ratio  is  slightly  greater  than  one  is  because 
the  transmission  coefficients  for  verticaJ  polarization  (426)  and  (456)  are  slightly 
greater  than  those  of  horizonteil  polarization  (42a)  and  (45a). 


^hh  (dB) 

‘’hvl'^hh 

P 

oblate 

4 

0,0625 

-24.6 

1.002 

0.108 

0.784 

oblate 

2 

0.25 

-24.2 

1.0006 

0.086 

0.829 

sphere 

1 

1 

-23.5 

1.0 

0.0 

1.0 

prolate 

Q.5 

4 

-25.1 

1.0001 

0.214 

0.572 

prolate 

0.25 

16 

-26.0 

1.0002 

0.265 

0.470 

Table  1.  Covariance  matrix  elements  for  five  types  of  scatterers.  Fre¬ 
quency  =  1.1  GHz, =  40°,  /  =  0.005,  <*  =  (30.8  -h  tl.8)eo, 
and  €{,  =  «(). 


Sliapc 

Ip  (mm) 

Ij  (mm) 

<^hh  (dB) 

‘'vv/chh 

‘'kvl<’hh 

P 

oblate 

4 

0  0625 

-24.8 

1.006 

0.107 

0.787 

oblate 

2 

0.25 

-24.4 

1.002 

0.085 

0.830 

sphere 

1 

1 

•23.5 

1.0003 

0.0 

1.0 

prolate 

0.5 

4 

-25.3 

1.001 

0.212 

0.577 

prolate 

0.25 

16 

•26.1 

1.002 

0.262 

0.477 

Table  2.  Covariance  matrix  elements  for  five  types  of  scatterers.  Fre¬ 
quency  =  1.1  GHz,  =  40°,/  =  0.01,  tj  =  (30.8 -h  tl.8)  eo, 
and  <(,  =  «o- 

The  nonspherical  geometry,  however,  has  significant  effects  on  o'hv/^hh  P- 
The  imaginary  part  of  p  is  not  shown  in  the  table,  since  they  are  very  close  to  zero. 
It  can  be  seen  that  when  the  scatterers  become  less  spherical,  <^hvl^hh  increases 
from  zero  and  p  decreases  from  1.  When  the  eccentricity  of  the  scatterers  gets 
higher,  the  returns  for  HH  and  VV  will  be  dominated  more  and  more  by  scatterers 
with  different  orientations  (i.e.,  different  location).  As  a  result,  the  correlation 
between  HH  and  VV  decreases.  As  a  matter  of  fact,  for  completely  randomly 
oriented  dipoles  both  the  ratio  and  p  are  equal  to  1/3.  This  is  consistent 

with  the  trend  shown  in  Table  1. 

The  results  shown  in  Tables  1  and  2  are  numerically  integrated  from  (53)  and 
(54).  Because  the  integrands  are  relatively  smooth  analytic  functions,  Simpson’s 
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rule  is  used  for  numericail  integration.  The  total  computation  time  for  one  co- 
variance  matrix  is  about  four  seconds  on  our  VAXstation  3500.  The  approximate 
formula  given  in  (57)  are  also  evaluated  and  the  approximate  results  are  within 
2%  error  of  the  exact  results  presented. 

The  parameters  <Tiiv/‘^hh  P  essentially  unchanged  when  the  fractional 
volume  of  the  scattercrs  changes  from  /  =  0.5%  (Table  1)  to  /  =  1%  (Table  2). 
Results  for  /  up  to  a  few  percent  have  aJso  been  evaluated  and  are  also  within 
a  few  percent  difference  with  the  numbers  shown  in  Table  1  for  those  two  pa¬ 
rameters.  This  indicates  that  (fhvl^hh  P  serve  as  useful  parameters  in 
identifying  and  distinguishing  the  geometry  of  scattercrs  for  the  remote  sensing 
of  vegetation  where  /  is  usually  smaller  than  a  few  percent. 

APPENDIX  A:  EFFECTIVE  PERMITTIVITY  UNDER  THE  LOW  FREQUENCY 
APPROXIMATION 

Assuming 

|(!e//)ol  <  1 

then,  from  (16), 

V/ -  +  {A2) 

Substituting  (A2)  and  (19)  into  (15),  then 

{D{f))  =  tgl(E(T))  +  j  df‘  (|(r)  .  PS^s(r,r^)  .  {(r^))  •  {W))  (A3) 

Assuming  the  frequency  is  low  enough  so  that  the  variation  of  mean  field  (P(f')) 
is  negligible  within  a  few  correlation  lengths,  we  can  approximate  (A3)  as 

(D(f))  =  +  klto  J  df  (!(r)  •  P$dg{f,f') .  !(r')))  •  (£(?))  (A4) 

Therefore,  we  can  define  the  effective  permittivity  of  the  medium  under  the  low 
frequency  approximation  as 

—  — 

«e//  -  *9-^  + ‘Ofe//  (A5) 

where 

th  =  kljd^  {1(f)  .  P5^j(f,f')  •!(f^))  (A6) 


APPENDIX  B:  EXPRESSION  FOR  ,  5.,  and  5, 

For  the  correlation  function  (37)  in  the  local  coordinate,  the  integration  of  a  more 
general  case  than  (27)  has  been  carried  out  analytically  by  Nghiem  [8].  In  this 
paper,  results  are  presented  for  isotropic  background  medium  only  and  are  given 
by 


1002 


Yueh  et  al. 


where 


ip 

^  =  u 

iB2a) 

(B2b) 

(52c) 

-1  =  kl  ^9  ‘1 

(B2d) 

Note  that  the  branch  cuts  of  In(-)  and  v/"  ^re  chosen  along  the  negative  real 
2Lxis  of  the  argument  in  the  above  expressions.  Taking  the  zero  frequency  limit  of 
(Bl),  we  obtain 


where 


Sp  =  ,  (l-fa^)tan  ^a  —  a 

^  2ega^ 

(B3a) 

<o(l  +  “^)  ^  -1 

Sz  = - ? -  a  —  tan  a 

{B3b) 

(54) 

APPENDIX  C:  EFFECTIVE  PERMITTIVITY  FOR  DISCRETE  SCATTERER  TH¬ 
EORY 

Under  the  low  frequency  assumption  and  in  the  limit  of  small  fractional  volume  for 
scatterers,  the  effective  permittivity  for  randomly  oriented  spheroids  for  discrete 
scatterer  theory  can  be  obtained  from  Equation  (10)  (pp.  534)  in  [16]  by  taking 
the  small  fractional  volume  limit.  The  result  for  spheroidal  scatterers  is  given  as 

^^1  =  1  + f  {2ya  +  yc)  +  /  (2yl  +  i  (Cl) 

where 

k  =  koy/^  (C2a) 


(C26) 

(C2c) 

(C3) 
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For  prolate  spheroids  (c  >  a) 


Ac  =  — 


1 


"  “  c3e3 


2e  +  In 


1  —  e 


1  +  e 


with 


e  =  A/l - =■ 


and  Aa  can  be  obtained  from  the  following  formula, 

2 


2i4fl  4"  Ac  — 


2 

a^c 


((74) 

((75) 

(C6) 
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1.1  Introduction 

With  advances  in  polarimetric  radar  technology,  polarimetry  has 
become  important  to  the  remote  sensing  of  geophysical  media.  Fully 
polarimetric  radar  signals  convey  additional  information  regarding  the 
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Figure  1.2.2  Geometrical  representatlona  of  polarisations. 
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